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Abstract 

We present a Hamiltonian framework for higher-dimensional vortex filaments (or membranes) and 
vortex sheets as singular 2-forms with support of codimensions 2 and 1, respectively, i.e. singular 
elements of the dual to the Lie algebra of divergence-free vector fields. It turns out that the localized 
induction approximation (LIA) of the hydrodynamical Euler equation describes the skew-mean- 
curvature flow on vortex membranes of codimension 2 in any R™ , which generalizes to any dimension 
the classical binormal, or vortex filament, equation in R 3 . 

This framework also allows one to define the symplectic structures on the spaces of vortex sheets, 
which interpolate between the corresponding structures on vortex filaments and smooth vorticities. 



Contents 

1 The vorticity form of the Euler equation 3 

1.1 The Euler equation on vorticity 3 

1.2 Smooth and singular vorticities [4 

2 Singular vorticities in codimension 2: point vortices and filaments 6 

2.1 Point vortices in 2D 6 

2.2 Vortex filaments in 3D |7 

3 Vortex membranes and skew- mean- curvature flow 

4 The localized induction approximation (LIA) in higher dimensions [ll| 

5 Singular vorticities in codimension 1: vortex sheets 

5.1 Vortex sheets as exact 2-forms 13 

5.2 Definition and properties of the symplectic structure on vortex sheets \\A 

6 Appendix: Derivation of the LIA in higher dimensions 

6.1 Generalized and localized Biot-Savart formulas 16 

6.2 Regularization of velocity 17 

6.3 Regularization of energy Il8 



'Department of Mathematics, University of Toronto, Toronto, ON M5S 2E4, Canada; e-mail: 
khesinOmath. tor onto . edu 



1 



Preface 



Vladimir Arnold's 1966 seminal paper [lj in which he introduced numerous geometric ideas into hy- 
drodynamics influenced the field far beyond its original scope. One of Arnold's remarkable and, in my 
opinion, very unexpected insights was to regard the fluid vorticity field (or the vorticity 2-form) as an 
element of the dual to the Lie algebra of the fluid velocities, i.e., the algebra of divergence-free vector 
fields on the flow domain. 

In this paper, after a review of the concept of isovorticed fields, which was crucial, e.g., in Arnold's 
stability criterion in fluid dynamics, we present an "avatar" of this concept, providing a natural frame- 
work for the formalism of vortex membranes and vortex sheets. In particular, we present the equation 
of localized induction approximation, which turns out to be the skew-mean-curvature flow in any di- 
mension. We also show that the space of vortex sheets has a natural symplectic structure and occupies 
an intermediate position between vortex filaments in 3D (or point vortices in 2D) equipped with the 
Marsden-Weinstein symplectic structure on the one hand and smooth vorticity fields with the Lie- 
Poisson structure on them on the other hand. 

Before launching into hydro dynamical formalism in this memorial paper, I would like to recall an 
episode with Vladimir Igorevich, related not to fluid dynamics, but rather to his equally surprising 
insights in real life: his remarks were always witty, to the point, and often mischievousQ 

Back in 1986 Arnold became a corresponding member of the Soviet Academy of Sciences. This 
was the time of "glasnost" and "acceleration" : novels of many formerly forbidden authors appeared in 
print for the first time. Jacques Chirac, Prime Minister of France at the time, visited Moscow and gave 
a speech in front of the Soviet Academy in the Spring of 1986. The speech was typeset beforehand 
and distributed to the Academy members. Arnold was meeting us, a group of his students, right 
after Chirac's speech and brought us that printout. Chirac, who knows Russian, mentioned almost 
every disgraced poet or writer of the Soviet Russia in his speech: it contained citations from Gumilev, 
Akhmatova, Mandelshtam, Pasternak... And on the top of this printout, above the speech, was the 
following epigraph in Arnold's unmistakable handwriting: 

"... H donyvny: ycnexu uatuu 6ucmpu 
Ho zde mc y uac Munucmp-deMazoz? 
IJycmb npo6epym ece chucku u pezucmpu, 
H numb py6jieu 6yM.aoK.Hux da,M e 3ajioz; 
Burnt Mootcem, ux eo Opawwuu neMajio, 
Ho ua Pycu ux nem - a He 6ueajio!" 

A.K. ToAcmou "Coh Honoea" (1873^ 
Who could take this speech seriously after such a tongue-in-cheek epigraph? As a curious aftermath, 
Arnold and Chirac shared the Russia State Prize in 2007. 

x At times it was hard to tell whether he was being serious or joking. For instance, when inviting the seminar participants 
for an annual ski trip outside of Moscow, Arnold would say: "This time we are not planning too much, only about 60km. 
Those who doubt they could make it — need not worry: the trail is so conveniently designed that one can return from 
interim bus stops on the way, which we'll be passing by every 20km." 

2 English translation by A.B. Givental: 

"... Our nation's rise is, I concur, gigantic, 
But demagogs among our statesmen?! Let 
all rosters, in the fashion most pedantic, 
be searched, I'd put five rubles for a bet: 
There could be more than few in France or Prussia, 
but are - and have been - none in mother Russia!" 
from "Popov's Dream" by A.K. Tolstoy (1873) 
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Returning to mathematics, I think that the ideas introduced by Arnold in pQ, so natural in retro- 
spect, are in fact most surprising given the state of the art in hydrodynamics of the mid-60s both for 
their deep insight into the nature of fluids and their geometric elegance and simplicity. In the next 
section we begin with a brief survey of the use of vorticity in a few hydrodynamical applications and 
discuss how it helps in understanding the properties of vortex filaments and vortex sheets. 



1 The vorticity form of the Euler equation 

Consider the Euler equation for an inviscid incompressible fluid filling a Riemannian manifold M 
(possibly with boundary). The fluid motion is described as an evolution of the fluid velocity field v in 
M which is governed by the classical Euler equation: 

d t v + (v, V> = -Vp. (1) 

Here the field v is assumed to be divergence-free (divv = 0) with respect to the Riemannian volume 
form n and tangent to the boundary of M . The pressure function p is defined uniquely modulo an 
additive constant by these restrictions on the velocity v. The term (v, V)i> stands for the Riemannian 
covariant derivative V v v of the field v in the direction of itself. 



1.1 The Euler equation on vorticity 

The vorticity (or Helmholtz) form of the Euler equation is 

d t £ + L v t; = 0, (2) 

where L v is the Lie derivative along the field v and which means that the vorticity field £ := curl v is 
transported by (or "frozen into" ) the fluid flow. In 3D the vorticity field £ can be thought of as a vector 
field, while in 2D it is a scalar vorticity function. In the standard 2D-space with coordinates (xi,X2) 
the vorticity function is curl v := dv2/dx\ — dvi/dx2, which can be viewed as the vertical coordinate 
of the vorticity vector field for the 2D plane-parallel flow in 3D. The fact that the vorticity is "frozen 
into" the flow allows one to define various invariants of the hydrodynamical Euler equation, e.g., the 
conservation of helicity in 3D and enstrophies in 2D. 

The Euler equation has the following Hamiltonian formulation. For an n-dimensional Riemannian 
manifold M with a volume form /j, consider the Lie group G = Diff^(M) of volume-preserving diffeo- 
morphisms of M. The corresponding Lie algebra g = Vect M (M) consists of smooth divergence- free 
vector fields in M tangent to the boundary dM: 

Vect M (M) = {V € Vect(M) | L vi x = QandV\\dM} E 

The natural "regular dual" space for this Lie algebra is the space of cosets of smooth 1-forms on M 
modulo exact 1-forms, g* = Q 1 (M) / dQ° (M) , see e.g. [HE]. The pairing between cosets [rj] of 1-forms 
rj and vector fields W G Vect^(M) is given by 

(M,W) := I iwVV, (3) 
Jm 

3 We usually denote generic elements of Vect fl (A/), as well as variation fields, by capital letters, while keeping the small 
v notation for velocity fields related to the dynamics. 



3 



where iw is the contraction of a differential form with a vector field W. The Euler equation ([T]) on the 
dual space assumes the form 

dt[Tj\+L v [rj\ = 0, 

where [rj] E Q 1 (M)/dfl°(M) stands for the coset of the 1-form r/ = related to the velocity vector 
field v by means of the Riemannian metric on M. (For a manifold M equipped with a Riemannian 
metric (., .) one defines the 1-form v" as the pointwise inner product with vectors of the velocity field 
v: V \W) := (v,W) for all W G T X M, see details in [Hi].) 

Instead of dealing with cosets of 1-forms, it is often more convenient to pass to their differentials. 
The vorticity 2-form £ := dv^ is the differential of the 1-form rj = v . Note that in 3D the vorticity 
vector field curl v is defined by the 2-form £ via i C ur\vH = £ for the volume form /i. In 2D curl v is 
the function curl v := The definition of vorticity £ as an exact 2-form in M makes sense for any 
dimension of the manifold M . This point of view can be traced back to the original papers by Arnold, 
see e.g. [21 [3] 

Such a definition immediately implies that: 

i) the vorticity 2-form £ := dij is well-defined for cosets [rj]: 1-forms r] in the same coset have equal 
vorticities, and 

ii) the Euler equation in the form ([2]) or dt(drj) + L v {drj) = means that the vorticity 2-form 
£ = drj is transported by (or frozen into) the fluid flow in any dimension. The latter allows one 
to define generalized enstrophies for all even-dimensional flows and helicity-type integrals for all odd- 
dimensional ideal fluid flows, which turn out to be first integrals of the corresponding higher-dimensional 
Euler equation, see e.g. [I]. This geometric setting can be rigorously developed within the Sobolev 
framework for H s diffeomorphisms and vector fields on M for sufficiently large s, see [6]. To present 
the geometric ideas and include singular vorticities we keep things formal in what follows. 

Remark 1.1. This point of view on vorticity was the basis for Arnold's stability criterion. Namely, 
steady fluid flows are critical points of the restriction of the Hamiltonian (which is the kinetic energy 
function defined on the dual space) to the spaces of isovorticed fields, i.e. sets of fields with diffeomorphic 
vorticities. If the restriction of the Hamiltonian functional has a sign-definite (positive or negative) 
second variation at the critical point, the corresponding steady flow is Lyapunov stable. This is famous 
Arnold's stability test. In particular, he proved (see e.g. [U Hj) that shear flows in an annulus with 
no inflection points in the velocity profile are Lyapunov stable, thus generalizing the Rayleigh stability 
condition. 



1.2 Smooth and singular vorticities 

Another consequence of such a point of view on vorticity, which is of the main interest to us, is the 
existence of the Poisson structure. 

Let M be an n-dimensional Riemannian manifold with a volume form /x and filled with an in- 
compressible fluid. As we discussed above, the vorticity of a fluid motion geometrically is the 2-form 
defined by £ := dv , where v is the 1-form obtained from the vector field v by the metric lifting of 
indices. Assume that i7 1 (M) = to simplify the reasoning below. Then the space of vorticities {£}, 
i.e. the space of exact 2-forms dil 1 (M), coincides with the dual space to the Lie algebra Vect At (M) of 
divergence-free vector fields. Indeed, Vect A1 (M)* ~ SV/dO ~ dQ 1 where the latter identification holds 
since iJ x (M) = 0. 
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Remark 1.2. As the dual space to a Lie algebra, the space of vorticities Vect^(M)* = {£} has the 
natural Lie-Poisson structure. Its symplectic leaves are coadjoint orbits of the corresponding group 
Diffw(M). Here such orbits are sets of fields with diffeomorphic vorticities on M, with the group action 
being the action of volume-preserving diffeomorphisms on vorticity 2-forms. The Euler equation defines 
a Hamiltonian evolution on these orbits. 

The corresponding (Kirillov-Kostant) symplectic structure on orbits in Vect /J (M)* is given by the 
following formula. Let V and W be two divergence-free vector fields in M, which we regard as a pair of 
variations of the point £ in Vect M (M)*. The Kirillov-Kostant symplectic structure on coadjoint orbits 
associates to a pair of such variations tangent to the coadjoint orbit of the vorticity £ the following 
quantity: 

uf K (V,W):=(d~H, [V,W]) = (n, [V,W]) = f r } M [V)W] ^= f £AW^- (4) 

Jm J m 

Here the 1-form rj = d _1 £ is a primitive of the vorticity 2-form £, and [V, W] is the commutator of the 
vector fields V and W in M. Note that for divergence-free vector fields V and W their commutator 
satisfies the identity iww}! 1 = d{iyiw^)-, which implies the last equality in (j3|). 

In this paper we deal with singular vorticities. Regular vorticities have support of full dimension, 
i.e. of codimension in M, while singular ones have support of codim > 1. Singular vorticities form 
a subspace in (a completion of) the dual space Vect M (M)* = d$7 1 (M). Note that since vorticity is a 
(possibly singular) 2-form (more precisely, a current of degree 2), its support has to be of codim < 2. 
(E.g., if support is of codim = 3, it corresponds to a singular 3-form. We refer to [8] for details on 
currents.) 

The most interesting cases of support are of codimension 1 (vortex sheets) and codimension 2 (point 
vortices in 2D, vortex filaments in 3D and vortex membranes, or higher filaments, for any dimension). In 
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the next sections we start with the codimension 2 case, and deal with the codimension 1 case towards 
the end of the paper. The main types of singular vorticities, as well as related to them symplectic 
structures and Hamiltonian equations studied below, are summarized in the table above. While the 
goal of this paper is partially expository, and various facts on vortex filament dynamics are scattered 
in the extensive literature, certain results presented below (in particular, the Hamiltonian framework 
for vortex sheets, the skew-mean-curvature flows and the LIA in any dimension) are apparently new. 



2 Singular vorticities in codimension 2: point vortices and filaments 
2.1 Point vortices in 2D 

Let M be the 2-dimensional Euclidean plane M 2 . Let the 2D vorticity £ be supported on N point 
vortices: £ = Ylf=i K j = Ylf=i K j ~ z j)' where Zj = (xj,yj) are coordinates of the jth point 
vortex in M 2 = C 1 with the standard area form \i = dx A dy. Kirchhoff's theorem states that the 
evolution of vortices according to the Euler equation is described by the system 

d U ■ d U 

W = — , KjV> = - Wj , 1<J<N. (5) 

This is a Hamiltonian system in M. 2N with the Hamiltonian function 

1 N 

j<k 

and the Poisson structure is given by the bracket 

if } = V— (— — -— — ^ (6) 

' !J fri K J \dxjd yj dyjdxjj ' 

One can derive the above Hamiltonian dynamics from the 2D Euler equation in the vorticity form 
9t(, = {ip, £,} j where £ is a vorticity function in M? and the stream function (or Hamiltonian) ip of the 
flow satisfies = £, see e.g. [10]. While the system (J5J) goes back to Kirchhoff, its properties for 
various numbers of point vortices and versions for different manifolds have been of constant interest, 
see e.g. [114 I10| . The cases of TV = 2 and N = 3 point vortices are integrable, while those of N > 4 
are not. The subtle issue of in what sense the equation of N point vortices approximates the 2D Euler 
equation as N — > oo is treated, e.g., in [T5]FI 

The origin of the Poisson bracket (jSJ) is explained by the following 

Proposition 2.1. (|14j) The Poisson bracket is defined by the Kirillov-Kostant symplectic structure 
on the coadjoint orbit of the (singular) vorticity £ = J^jLi K j ~ z j) * n completion of) the dual 
of the Lie algebra q = Vect^(R 2 ) of divergence free vector fields in M 2 . 



4 The Hamiltonian system for the point vortex approximation of the 2D Euler equation is reminiscent of the Calogero- 
Moser system for the evolution of poles of rational solutions of the KdV equation. 
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Proof. Indeed, a vector tangent to the coadjoint orbit of such a singular vorticity v can be regarded 
as a collection of vectors {Vj} in M. 2 attached at points Zj. Then for a pair of tangent vectors the 
corresponding Kirillov-Kostant symplectic structure becomes the weighted sum of the corresponding 
contributions at each point vortex Zj with strengths Kj as the corresponding weights: 

uf K (V, W) := f £ A i v i wl i = V Kj n{V s , Wj) . 

The Poisson bracket ©, being the inverse of the symplectic structure, has the reciprocals of the weights 
Kj . □ 



2.2 Vortex filaments in 3D 

By passing from 2D to 3D we move from point vortices to filaments. Vortex filaments are curves in R 3 
being supports of singular vorticity fields. They are governed by the Euler equation 

$e + A,£ = o, (7) 

where v = curl _1 £ and the vorticity field (or a 2-form) £ has support on a curve 7 CK 3 . (Note that the 
exactness of the form £ implies that 7 is a boundary of a 2-dimensional domain, i.e., in particular, its 
components are either closed or go to infinity.) The Euler dynamics is nonlocal in terms of the vorticity 
field, or 2-form, £ since it requires finding the field v = curl £. 

The localized induction approximation (LIA) of the vorticity motion is a procedure which allows 
one to keep only the local terms in the vorticity Euler equation, as we discuss below. In M 3 the 
corresponding evolution is described by the vortex filament equation 

da = i x 7 " , (8) 

where j(-,t) C M 3 is a time-dependent arc-length parametrized space curve. For an arbitrary parametriza- 
tion the filament equation becomes £^7 = k ■ b, where k and b = t x n stand, respectively, for the 
curvature value and binormal unit vector of the curve 7 at the corresponding point. This equation is 
often called the binormal equation. (The equivalence of the two equation forms is straightforward: for 
an arc- length parametrization the tangent vectors t = dj/88 = 7' have unit length and the acceleration 
vectors are 7" = dt/d9 = k ■ n, i.e. dtj = 7' x 7" becomes dtj = k • b, where the latter equation is 
valid for an arbitrary parametrization.) 

Remark 2.2. Here we briefly recall the LIA derivation of the equation (JSJ) in 3D, see e.g. {5], and give 
more details in Section H] and Appendix. 

Assume that the velocity distribution v in ]R 3 has vorticity £ = curl-u concentrated on a smooth 
embedded arc- length parametrized curve 7 C I 3 of length L. Then 

=c £ 5(q- 7 (9,t))^de. 

Here 5 is the delta- function in R 3 and the constant C, the strength of the filament, is the flux of £ 
across (or, which is the same, the circulation of v over) a small contour around the core of the vortex 
filament 7. Note that the exactness of the 2-form £ also implies that the filament strength C is indeed 
constant along 7. 
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The Biot-Savart law allows one to represent the velocity field in terms of its vorticity: 



4*Jm \k-i\\ 3 ^All«-7(«,0II 3 *> 

By utilizing the fact that the time evolution of the curve 7 is given by the velocity field v itself: 
^(0,t) = v( r f(6,t),t) we come to the following integral: 

fr (M ) = -£/ ->(»■ 0-7(j.O x g7^ 

» V 4iry 7 || 7 (e,t)-7(«,0ll 3 9« 

This integral is divergent with the main singularity coming from the points on the curve 7 close to each 
other on the curve (i.e. with small \0 — 0\). Given the Taylor expansion of j(0) in (0 — 0) yields 



&1 (0 t)-° 



80 80* 1 <" "< +U[1) 







Since the right-hand side is divergent, for a small e we take the truncation of the integral by considering 
only the part \0 — 0\ > e of the integration domain [0, L]. The corresponding integral is of order hie. 
Rescaling time by means of t — > i(C/87r)lne we are keeping only the singularity term and neglecting 
others as e — > 0. 

This way we obtain the vortex filament equation ([8]). It is also called the localized induction 
approximation (LIA) since the velocity field d^y/dt of the curve 7 is induced by its own vorticity, i.e. 
vorticity supported on the curve, while only parts of the curve sufficiently close to a given point 7(0) 
determine the velocity field at that point. (We discuss the above limit in higher dimensions in the next 
two sections.) Note that in 2D point vortices interact with each other but not with themselves (as 
manifested by the Kirchhoff Hamiltonian) , i.e. the localization in 2D would give the zero LIA equation. 

Remark 2.3. This binormal equation is known to be Hamiltonian relative to the Marsden-Weinstein 
symplectic structure on non-parametrized space curves in K 3 . Recall that the Marsden-Weinstein 
symplectic structure is defined on oriented curves 7 by 



UJ. 
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(V,W):= / %y% wi i= \ n(V,W,j')d0 (10) 



where V and W are two vector fields attached to the curve 7 and regarded as variations of this curve, 
while the volume form [i is evaluated on the three vectors V, W and 7' = 8^/80. One can see that this 
integral does not depend on the parameter on the curve j(0). 

Equivalently, this symplectic structure can be defined by means of the operator J of almost complex 
structure on curves: any variation, i.e. vector field attached at the oriented curve 7, is rotated by the 
operator J in the planes orthogonal to 7 by ir/2 in the positive direction (which makes a skew-gradient 
from a gradient field) , see details in |14|, S] . 

One can show that this is the Kirillov-Kostant symplectic structure on the coadjoint orbit of the 
vorticity £ 7 supported on the curve 7 and understood as a point in a completion of the dual of the 
Lie algebra: £ 7 € cZQ 1 (M 3 ) = Vect^(R 3 )*. The pairing of 7 and a divergence- free vector field V can be 
defined directly as (7, V) := Flux V| CT , where a is an oriented surface whose boundary is 7 = da. 

Remark 2.4. As discussed above, the Euler equation ([7]) is Hamiltonian with the Hamiltonian function 
given by the kinetic energy. The energy E(v) = \ f M (v,v) fi is local in terms of velocity fields, but it 
is nonlocal in terms of vorticities: E(£) = | f M (curl -1 £, curl^ 1 ^) [i. It turns out that after taking the 
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localized induction approximation, when we keep only the local terms, the filament equation remains 
Hamiltonian with respect to the same Marsden-Weinstein symplectic structure, but with a different 
Hamiltonian (see Sections EJ H] and Appendix) . 

The corresponding new Hamiltonian functional turns out to be the length functional of the curve: 
#(7) = length(7) = J \\ r y'(6)\\d6, see e.g. [4]. Indeed, the variational derivative, i.e. the "gradient," 
of this length functional H is 5H/5^f = —7" = — t' = —k ■ n, where t and n are, respectively, the unit 
tangent and normal fields to the curve 7. The dynamics is given by the corresponding skew-gradient, 
which is obtained from 5H/5^ by applying J for the above symplectic structure. This operator, rotating 
the plane orthogonal to t by vr/2, sends — k ■ n to k ■ b. In the next section we discuss how the Marsden- 
Weinstein symplectic structure and binormal equation are generalized to higher dimensions. 

The LIA evolution is close to the actual Euler evolution of a vortex filament only for a short time 
(when the local term in dominant). For large times the LIA filament may, e.g., self-intersect, while the 
incompressible Euler dynamics has a frozen-in vorticity and it does not allow topology changes of the 
filaments. 



3 Vortex membranes and skew-mean-curvature flow 

For a smooth hypersurface in the Euclidean space M n its mean curvature is a function on this surface. 
Similarly, one can define a mean curvature vector field for a smooth submanifold of any dimension I. 

Definition 3.1. a) Let P be a smooth submanifold of dimension I in the Euclidean space M n . Its 
second fundamental form at a point p G P is a map from the tangent space T p P to the normal space 
N p P. The mean curvature vector MC(p) G N p P is the normalized trace of the second fundamental 
form at p, i.e. the trace divided by /. 

b) Equivalently, the mean curvature vector MC(p) G N p P is the mean value of the curvature vectors 
of geodesies in P passing through the point p when we average over the sphere S 1-1 of all possible unit 
tangent vectors in T p P for these geodesies. 

Now consider a closed oriented embedded submanifold (membrane) P of codimension 2 in R n 
(or more generally, in a Riemannian manifold M n ) with n > 3 and the Marsden-Weinstein (MW) 
symplectic structure on such submanifolds. Recall that the Marsden-Weinstein symplectic structure 
oj mw on membranes of codimension 2 in M n (or in any n-dimensional manifold) with a volume form // 
is defined similar to the 3-dimensional case: two variations of a membrane P are regarded as a pair of 
normal vector fields attached to the membrane P and the value of the symplectic structure on them is 

CO^ W (V,W) := j i V i W p. 

Here iyiwH is an (n — 2)-form integrated over P. Note that this symplectic structure can be thought of 
as the "total" averaging of the symplectic structures in each normal space N p P to P. (The Marsden- 
Weinstein structure in higher dimensions was studied in pQ[§]-) 

Now we define the Hamiltonian function on those membranes by taking their (n — 2)-volume: 
H(P) = volume(P) = J p /J-p, where fip is the volume form of the metric induced from M. n to P. For 
instance, for a closed curve 7 in M 3 this Hamiltonian is the length functional discussed in Remark 
12.41 Note that to define the MW structure one only needs the volume form on W 1 , while to define the 
Hamiltonian one does need a metric. 
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Theorem 3.2. In any dimension n > 3 the Hamiltonian vector field for the Hamiltonian H and the 
Marsden- Weinstein symplectic structure on codimension 2 membranes P C MJ 1 is 

v H (p) = C n -J(MC(p)), 

where C n is a constant, J is the operator of positive tt/2 rotation in every normal space N p P to P, and 
MC(p) is the mean curvature vector to P at the point p. 

This statement, as well as the proof below, is valid for any Riemannian manifold M. The expression 
of vh via the trace of the second fundamental form without reference to the mean curvature appeared 
in [9], Proposition 3. For 4D this theorem was obtained in [16] . Here and below we use the notation 
C n for some constant depending on the dimension in the case of M n , or on the geometry of M n in the 
general case, but not on the membrane P. In the theorem above C n = 4 — 2n. 

Proof. Since the MW symplectic structure is the averaging of the symplectic structures in all 2- 
dimensional normal planes N P P, the skew-gradient (i.e. the Hamiltonian vector) for any functional 
on submanifolds P is obtained from its gradient field attached at P by the application of the almost 
complex structure J. The latter is the positive rotation by ir/2 in each normal plane. (Orientations of 
M. n and P determine the orientation of N p P and hence the positive direction of rotation in N p P is well 
defined.) Thus to prove vh{p) = const • J(MC(p)) we need to show that the gradient, i.e. the first 
variation, of the volume functional H{P) is 

5H 

— (p) = const • MC(p) . 
oP 

On the other hand, the fact that the mean curvature vector field is the gradient for the volume 
functional is well-known, see e.g. [7|. A quick argument follows from the observation that for a 
variation Pt defined by a normal vector field W attached at P of dimension / the volume changes at 
the rate 

j t H(P t ) = -2l J (W,MC)hp. 

The latter equality can be verified for a variation confined to a local chart parametrizing a neighborhood 
of a point p € P. Let di be coordinate unit vectors in this chart and <f> is the chart parametrization 
map. Then the induced metric on P in local coordinates around the point p is gij = (<f>*di,<f)#dj) and 
the volume variation is d/dtdet(gij). By choosing the coordinates so that gijip) = Sij at t = one has 

d/dtdet(gij) = \x(d/dtgij) = tr{Lw9ij) = 2ti(Vw<P*di,(j)*dj) = 2tr(V < £„g i W / , . 

Then using integration by parts one has 

d 



dt 



H{P t ) =2 J tr(V^ Si PM*d,)MP = -2^tr(W 5 V^*d,')MP = -21 J(W, MC) (J, P , 



since MC := ti('V ( j )it g i (j) t ,dj)/l. By applying this to P of dimension n — 2 one obtains C n = —2(n — 2) = 
4 - 2n. □ 

The mean curvature flow is often used to construct minimal surfaces in M 3 . For hypersurfaces it is 
directed along the normal, its value is given by the mean curvature, and it minimizes the hypersurface 
volume in the fastest way. 

Definition 3.3. The higher vortex filament equation on submanifolds of codimension 2 in R n is given 
by the skew-mean- curvature flow: 

d t P(p) = - J(MC(p)) . (11) 
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Note that the skew-mean-curvature flow introduced this way differs by the 7r/2-rotation from the 
mean-curvature one. Respectively, it does not stretch the submanifold while moving its points orthog- 
onally to the mean curvatures. In particular, the volume of the submanifold P is preserved under this 
evolution, as it should, being the Hamiltonian function of the corresponding dynamics. 

Remark 3.4. For dimension n = 3 the mean curvature vector is the curvature vector k ■ n of a curve 7: 
MC = Ami, while the skew-mean-curvature flow becomes the binormal equation: dtj = —J(k-n) = k-h, 
which for arc- length parametrization is given by the equation dt'j = 7' x 7". Unlike the case n = 3, for 
larger n > 4 the skew-mean-curvature flow is apparently non-integrable. 

Problem 3.5. a) Is there an analogue of the Hasimoto transformation for any n relating the higher vor- 
tex filament equation with the higher-dimensional (and already non-integrable) nonlinear Schrodinger 
equation (NLS)? 

b) Is there an analogue for any n of the gas dynamics equation equivalent to the vortex filament 
one in 3D, see [I]. 

c) Are there integrable approximations of the higher filament equation (llip . 

For n = 4 the question a) was posed in [16 . Recall that for n = 3 at any time t the Hasimoto 
transformation sends a curve 7(0) with curvature k{9) and torsion t(9) to the wave function ip{9) = 
k{9) exp{i f t(C) d(} satisfying the 1-dimensional NLS: idtip + ip" + ^^ip = 0. 

4 The localized induction approximation (LIA) in higher dimensions 

Let P n ~ 2 c M n , n > 3 be a closed oriented submanifold of codimension 2. Consider the vorticity 2-form 
£p supported on this submanifold: £p = C -5p. We will call P a higher (-dimensional) vortex filament or 
membrane. Note that the exactness (or closedness) of the 2-form £p implies that the membrane strength 
C is constant, while the integrals of £p over 2-dimensional surfaces with boundary not intersecting P 
are well-defined and depend only on the homology class of the boundary in the complement to P. 

We would like to find the divergence-free vector field v which has a prescribed vorticity 2-form £, i.e. 
£ = dv^ £ fi 2 (M n ). In dimension 3, where vorticity can be regarded as a vector field, the corresponding 
vector potential v in R 3 is reconstructed by means of the Biot-Savart formula ([9]). Now we are looking 
for its analogue in any dimension n > 3. The statements we are discussing in this section were obtained 
for n = 4 in [16J. The extentions of proofs to the general case are presented in Appendix. 

The singular <5-type vorticity 2-form £p is completely defined by the submanifold P. Denote by 
G(q,p) the Green function of the Laplace operator in M n , i.e. given a point q € 1" one has A p G(q,p) = 
Sg(p), the delta- function supported at q. 

Theorem 4.1. For any dimension n > 3 the divergence-free vector field v in M. n satisfying curlv = £p 
is given by the following generalized Biot-Savart formula: for any point q P one has 



where Proi N 'VpG(-,p) is the orthogonal projection of the gradient V p G(-,p) of the Green function G(-,p) 
to the fiber N p P of the normal bundle to P at p G P, the operator J is the positive rotation around p 
by tt/2 in this 2-dimensional space N p P, and [ip is the induced Riemannian (n — 2)-volume form on 
the submanifold P C W 1 . 
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In other words, 

v(q) := C n • ^sgrad p (G(q,p)\ NpP ) /xp(p) , 

by using the symplectic structure in N p P. Here G(q,p)\N p p is the restriction of the function G(q,p) 
to the normal plane N p P. These formulas use the affine structure of M. n , since in the integral averages 
vectors over P and attaches the total at the point q. In the case of an arbitrary manifold M n the 
Biot-Savart formula is more complicated, and we will use a round-about way to obtain the LIA for any 
M n , see Remark 14.51 below. 



Note that as the point q approaches the membrane P the vector field v(q) may go to infinity. 
Consider the following truncation of the integral above. For q € P and given e > take the integral 
over P for all points p satisfying \\q — p\\ > e, i.e. at the distance at least e from q: 

v t (q) :=C n - J J (Proj^V p G(q,p)) fj, P (p) . 

p&P, \\q— 

Theorem 4.2. (cf. [16] for 4D) For any dimension n > 3 the velocity field v defined in Theorem 
has the following asymptotic of the truncation v t : for g£Pc M. n one has 

lim^M = C n .J(MC( 9 )) . 
e^o lne 

By reparametrizing the time variable t — >• — (C n • lne)t to absorb the logarithmic singularity we 
come to the following LIA equation for a higher filament P C M n . 

Corollary 4.3. The LIA approximation for a vortex membrane (or higher filament) P in W 1 coincides 
with the skew-mean- curvature flow: 

d t P(q) = -J(MC(q)) , 

where MC(g) is the mean curvature vector at q € P. In particular, the LIA equation is Hamiltonian 
with respect to the Marsden-Weinstein symplectic structure and Hamiltonian function given by the 
volume of the membrane P. 

Consider now the energy Hamiltonian E{v) = \ f M {v, v) for M = W n and fast decaying divergence- 
free velocity vector fields v. As before, let £ be the vorticity 2- form of the field v, i.e. £ = dv . If 
the vorticity £p is supported on a membrane P C W 1 of codimension 2, the corresponding energy 
E{y) = \ Jw n { v i v ) M f° r the velocity v defined by cuilv = £p is divergent and requires a regularization. 
Consider the regularized energy 

E e(v) := - I (v,v e )fi. 



Theorem 4.4. (cf. [16] for 4D) For any dimension n > 3 the regularized energy E e (v) for the velocity 
of a membrane P C M n has the following asymptotics: 



lim = C n ■ [ up = C n ■ volume (P) . 

e^o lne Jp 



We refer to Appendix and [16] for details on the proofs for R n . As we discuss in Appendix, this 
regularization is also valid for any Riemannian manifold M. 



12 



Remark 4.5. When one passes from smooth to singular vorticities supported on membranes of codi- 
mension 2 the Euler dynamics requires regularization. Correspondingly, so does the associated energy 
Hamiltonian. On the other hand, the corresponding symplectic structure on smooth vorticities natu- 
rally descends to the MW symplectic structure on submanifolds (this is how it was defined in [TJ]) and 
does not need a regularization. 

This consistency explains why the hydro dynamical Euler equation remains Hamiltonian under the 
localized induction approximation. Indeed, the LIA takes the Hamiltonian Euler equation into the 
Hamiltonian skew-mean-curvature equation by "keeping only the logarithmic divergences" given by the 
local terms. 

For any manifold M the above consistency can be taken as the definition of the regularized dynamics, 
defined in Theorem 14.21 and Corollary 14.31 Namely, one can employ only the MW symplectic structure 
and regularization of the Hamiltonian, which uses only local properties of the Green function that hold 
for any M, in order to find the vortex dynamics in the general case. 



5 Singular vorticities in codimension 1: vortex sheets 
5.1 Vortex sheets as exact 2-forms 

Now we return to an arbitrary manifold M (with H X {M) = 0), but consider singular vorticities sup- 
ported in codimension 1. Introduce the following 

Definition 5.1. Vortex sheets are singular exact 2-forms, i.e. 2-currents of type £ = a A <5r, where 
r™ _1 C M n is a closed oriented hypersurface in M, Sr is the corresponding Dirac 1-current supported 
on r, and a is a closed 1-form on T. 

For a singular 2-form £ = a A 5-p to be exact either 

i) the closed 1-form a must be exact on T, i.e. a = df for a function / on T, or 

ii) the hypersurface V must be a boundary of some domain d~ l T C M and the closed 1-form a has 
to admit an extension to a closed 1-form a on d V. (For instance, T is a torus in M 3 while a = dO 
with 9 being one of the generating angles of the torus.) Note that under the assumption H 1 (M) = a 
closed hypersurface T is always a boundary. 

The above options come from the interpretation of £ = a A 5r- one can choose either a or 5r to be 
exact, while the other form closed, for the wedge product to be exact. 

For an exact form a = df the vortex sheet is fibered by levels of the function /. If a is a closed 
1-form, it is a function differential only locally, and the integral submanifolds of ker a foliate T. Thus 
the vortex sheets are fibered into filaments (of codimension 1 in T) in the former case and foliated in 
the latter. 

Example 5.2. If a is supported on a single hypersurface 7 in T (i.e. on a curve 7 C T for n = 3), then 
the vortex sheet £ = a A 5r = <5 7 reduces to the vorticity of the filament 7 C T. 

Remark 5.3. The corresponding primitive 1-forms ry satisfying £ = dr/ for the singular vorticity 2-form 
£ = a A 8r follows. 

i) For an exact a = df take 77 = f5r- 

ii) For a closed 1-form a extendable to a closed 1-form a on a domain <9 -1 T take as a primitive 
77 = the 1-form r\ = —Xd^T ' <5> where d~ l T is a domain bounded by the hypersurface T and Xd^T 
is its characteristic function. Indeed, 

dr] = —d(xd~ 1 r ' ot) = ~dXd~ 1 r A a = — 5r Aa = a/\5r=(,- 
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Note that the 1-form a and the domain d T are not defined uniquely, and this ambiguity corre- 
sponds to the ambiguity in the definition of a primitive 1-form n = d~ 1 !^. 

Vortex sheets £ understood as singular currents can be regarded as elements of a completion of the 
dual space Vect At (M)*. (It is convenient to change the order in this wedge product to £ = 5r A a in order 
to avoid the signs depending on the dimension of M in the pairing and symplectic structure below.) 

Definition— Proposition 5.4. The pairing of vortex sheets (i.e. singular vorticity currents) £ = 5^ A a 
with vector fields V £ Vect^(M) is defined by (cf. the pairing (EpJ 



where d 1 (5r A a) is a primitive 1-form for the vorticity £ = <5r A a. The pairing is well defined, i.e. it 
does not depend on the choice of d^ 1 . 

Proof. Indeed, 



Since H l (M) = the closed (n — l)-form iy[i is exact, and its primitives d may differ by an 

exact (n — 2)-form Then the form a A £ is exact on T and the corresponding pairing difference given 



For instance, for an exact a = df the pairing reduces to (d 1 (5r A df), V) = Flux (/V)|r- 

Remark 5.5. Suppose that T is the oriented boundary between two different parts Mj with velocity 
fields vi,V2 that are divergence- free and vorticity- free (i.e. locally potential flows). The vorticity is 
infinite at the interface T and here we describe how to define the 1-form a in the corresponding vortex 
sheet £ = a A Sr ■ 

Given a Riemannian metric on M, we prepare the 1-form Vj on Mj corresponding to the velocity 
Vj, respectively. Note that the forms in must be locally exact, m = dhj since curlfj = on Mj. 
Then locally a := {dh\ — dh,2)\r = dfi ~ df2- One can also define this 1-form a = d(fi — fe) by 
means of the vector field vr inside this vortex sheet T by using the metric restricted to T: locally 
vr := (d(f\ — /b))' = Proj|r(vi — ^2)- The proper sign of vr or the form a depends on the orientation of 
T: the latter defines the orientation of the corresponding exterior normal and hence signs of the fields 
v\ and V2 in this difference. 



5.2 Definition and properties of the symplectic structure on vortex sheets 

There is a natural symplectic structure on vortex sheets coming from the Lie-Poisson structure on 
Vect^(M)*. It extends the Marsden-Weinstein symplectic structure for filaments in M 3 and for mem- 
branes of codimension 2 in M n . The corresponding symplectic leaves are defined by isovorticed fields, 
i.e. fields with diffeomorphic singular vorticities a A 5-p. The corresponding symplectic structure on 
spaces of diffeomorphic vortex sheets is defined as follows. 





by the integral over Y is zero. 



□ 
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Definition 5.6. Given two vector fields V, W attached at T define the symplectic structure on variations 
of vortex sheets £ = 5r A a, i.e. pairs (T,a), by 



us r /\a{V,W) := J aAiyiw^- 

Theorem 5.7. The form cos r Aa coincides with the Kirillov-Kostant symplectic structure uj^ k on the 
coadjoint orbit containing the vortex sheet £ in Vect*(M). 

Proof. Adapt the formula Q for Kirillov-Kostant symplectic structure on the coadjoint orbit of £ to 
the case of a vortex sheet £ = 5r A a. Let V and W be two variations of £ given by divergence-free 
vector fields on M. Then by using the identity iww]^ = d{iyiw^) valid for divergence-free fields and 
specifying to the case of £ = 5r A a one obtains 



U) 



KK i 



(V, W) := / d x £ A i[v t w]V = / d x £ A d(iyi w \i) = / £ A iyiwl 1 



M JM JM 

= I 8r A a A iyiwH = / Oi A iyiwH = us rAa (V, W) . 
Jm Jr 



□ 



Remark 5.8. If a is supported on a curve 7 C T, i.e. £ = 5 7 , then co^(V,W) := J iyiw^- F° r 
a curve 7 C M 3 this is exactly the Marsden-Weinsten symplectic structure w on filaments, i.e. 
non-parametrized curves in M 3 , see (|10p . 



The evolution of vortex sheets £ = a A 5r is defined by the classical Euler equation in the vorticity 
form c\£ + L v £ = 0, where £ = curlu = cfoA This equation is Hamiltonian with respect to the above 
symplectic structure oj^ k . The standard energy Hamiltonian E(v) = ^ f M (v,v) fi defines a non-local 
evolution of the vortex sheet, similarly to the case of membranes. 

Let (/, 9) be coordinates on a vortex sheet a A 5r in R 3 where the exact 1-form a = df and the 
surface T is fibered into the filaments Tf being levels of the function /. The rough LIA procedure 
similar to the one described in Section [2] under the cut-off assumption e < \9 — 9\ < \f — f\ 2 leads 
to the binormal type equation: dtT = Tq x Tog , which is Hamiltonian with the Hamiltonian function 
H(T) := J length(r^) df. The latter may be understood as a continuous family of binormal equations. 
One may hope that other assumptions on the cut-off procedure lead to more interesting approximations. 



Problem 5.9. Describe possible analogues of localized induction approximations (LIAs) and the length 
Hamiltonian for vortex sheets. 

The Euler evolution of vortex sheets is described in the closed form by the Birkhoff-Rott equation, 
see e.g. |12j . The motion of vortex sheets is known to be subject to instabilities of Kelvin-Helmholtz 
type which lead to roll-up phenomena. It would be interesting to obtain this instability within the 
Hamiltonian framework for vortex sheets described above, cf. |13j . 



6 Appendix: Derivation of the LIA in higher dimensions 

In this Appendix we outline, following [TB] and extending it to any dimension, the generalized Biot- 
Savart formula and regularized energy for the vector fields whose vorticity is confined to membranes, 
i.e. submanifolds of codimension 2. 
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6.1 Generalized and localized Biot-Savart formulas 

Let v be a vector field in the Euclidean space W 1 , n > 3. Assume this field to be divergence free: 
divu = or, equivalently, d*v } = for the 1-form on R n . Its vorticity is the 2-form £ = dv . We are 
looking for a generalized Biot-Savart formula which would allow one to reconstruct the velocity field v 
for a given vorticity 2-form £, and in particular, for a given singular vorticity £p = 5p supported on a 
compact membrane P. 

Consider d*£ = d*dv^ = Av b = (Av) b . Then component-wise one has the Poisson equation on v: 
Avi = *(dxi A *d*£). 

Let G(-,p) be the Green function for the Laplace operator in W 1 B p. Then at any q £ M. n the 
components of the field-potential are 

Vi (q)= [ G(q,p)A*(d Xi A*d*0(p)Kp) = - [ (<9j, (*(d p G(q,p) A *£))") MO 5 ) > 

where di are the coordinate unit vectors in IR n . The vector field-potential v itself is 

v(q) = - [ (*(d p G(q, p) A M(P) , (12) 

JR™ 

which is the generalized Biot-Savart formula in the case of smooth vorticity £. 



Theorem 6.1. (= 14.11 ) For any dimension n the divergence- free vector field v satisfying curl-u = £p 
(i.e. dv^ = £p) is given by the following localized Biot-Savart formula: for any point q P one has 



j(q) := / sgrad p (G(q,p)\ NpP ) ^ P {p) 
J P 



(13) 



where G(q,p)\]y p is the restriction of the function G(q,p) to the normal plane N p P C M n . 



Proof. In order to set ^ to be = <5p we think of the latter in terms of local coordinates. Let 
ti,...,t n -2 be local coordinates along P, while v\,v% are coordinates normal to P near p £ P. Then 
locally £p(q) = C5 p (q) dv\ f\dv 2 where 5 p (q) is a delta-function supported at p € P, i.e. £p is the 5-type 
2-form in the transversal to P direction, and one has *£p = /ip. 

Then for q G M. n and p € P one has 

v(q) = ~ S p (q)(*(dpG(q,p) A fip)^ n(p) 

J R u V oui dv 2 J 

= ~ j P dV2 ~ ~^u^ dU ^) = f p s ^ ad p( G ( q >P)\ N P p } ^(P) ' 

where the last equality is due to the following (all derivatives of G(q,p) are in p, so we skip the index): 

dG dG \ ^ (9Cr _ dG - /^fi \ t /t-» ■ «— r ^-y\ 

□ 
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Remark 6.2. For q £" P the integrand expression above is smooth, since so is G(q,p) as a function 
of p € P. For q £ P the integral (|13p is well defined provided that the integration over P is replaced 
by that over P t = {p £ P | ||p — g|| > e}. As p — > q £ P the Green function has a singularity 
G(q,p) = C n |jr|j 2_n where r := p — q £ W 1 . Hence V P G = C n r/||r|| n , and therefore the integral is 
divergent. (Recall that C n stands for any constant depending on n.) This divergence is "local" in the 
sense that the contributions from p £ P close to q £ P make the velocity v(q) divergent, and this local 
contribution into v (q) is exactly what the LIA takes into account. 



6.2 Regularization of velocity 

Given e > consider a geodesic ball U t in the membrane P of radius e around a point q £ P. Define 
now a truncation v e (q) by integrating in f|13f) over P e := P\ U e instead of over P: 



Ve(q) ■= sgracL (G(q,p)\ NpP ) p, p . 



Theorem 6.3. (= 14.21 ) For any dimension n the approximation v e (q) has the following asymptotics: 
at any point q £ P one has 

Um = Cn . J (MC(?)) , 
e->o In e 

where MC(ij) is the mean curvature vector of the membrane P at q and the constant C n depends on n 
only. 

Proof. In order to find the asymptotics of how v t (q) — > oo as e — > we localize v e , i.e. confine the 
integration to a punctured neighborhood U e ^ a = {p £ P t | e < \\p — q\\ < a} C P e , since the integral 
outside of it, over P a = P e \ U t ^ a is finite. 

Set the origin of IR n at q, denote the radius vector from qtophyr = p — q £ M. n . Introduce 
the geodesic radial coordinate p and spherical multi-coordinate G in the ball U a of radius a inside the 
membrane P. Note that for the volume form on P of dimension n — 2 one has pp = p n ~^ dpd@. 

Then for a point p £ U e ^ a one has G(q,p) = C n /||r|| n_2 ~ C n /p n ~ 2 for the Green function, where ~ 
stands for the leading term in the corresponding expansion. Hence, V p G(q,p) ~ C n r/p n . Denote by 
u%, V2 normal coordinates to the codimension 2 membrane P near q. We have 

«e(9)~ / J(^oi N (V p G(q,p)))p P (p) = C n J pt^MllK^^^e. 

Now we fix (temporarily suppressing this notation) and denote by t{p) = dr/dp the tangent vector 
to the geodesic in direction 0: expand the following quantities in p near p = in the punctured 
neighborhood U t)Cl as follows: 

d^p)=d,M+P^)+0(p 2 ); 

r(p) = p|(0) + £0(0) + CV) = pt(0) + £|i(o) + 0(p 3 ) . 

Then for a given using (t, ^-) = — (0, c^J, which is implied by (t, d Ui )(p) = for any p, one obtains 
the following expansion of (r, d Vl )d Vl + (r,d U2 )d V2 in p: 



^ 2 



(0) + O(p 3 



17 



(O) + 0(/> 3 ) = -^curv t (O) + 0(p 3 ). 

Here curvt(O) is the vector of the geodesic curvature for the direction t at p = and fixed O, i.e. at 
the point £ R n , which stands for q. Restoring the dependence on we have 

vM ~cJ r j P 2 curv,(0,e ) p „_ 3dH@ 

JS»-3 Je 9 

= C n f ^ • J / curv t (O,0)d0 ~ C n -hie- J(MC(g)) 
by Definition 13.1b ) of the mean curvature vector. □ 



P_ 

' 2 



at 

dp 1 



, d Vl d, 



at 

dp' 



, d V2 d, 



'v 2 



6.3 Regularization of energy 

Obtain now a regularized expression for the corresponding energy of the velocity field v. Recall that 
the kinetic energy of a fluid moving with velocity v in a manifold M with a Riemannian volume form 
p is E(v) = \ f M (v, v)p = \ f M v* A *v b . 

Let £p be the vorticity 2-form supported on a membrane P C M n of codimension 2. As we will 
see below, the corresponding energy E(v) = \ f M (v,v) p for the velocity v satisfying curlv = £p is 
divergent. Following [16] define the regularized energy 

E e(v) = \ (v,V e )p. 

Theorem 6.4. (= 14.41 ) For any dimension n the regularized energy E e {v) has the following asymp- 
totics: 



lim = C n ■ [ pp = C n ■ volume (P) 

e^o me Jp 



Proof. First for any vector field v we rewrite the energy E(v) via vorticity by introducing the form- 
potential: 

E(y) = - I v* A *d* p = \ I £ A */3 
2 Jm 2 J M 

for the closed 2-form /3 satisfying d*f3 = v or, equivalently, A/3 = dd*(3 = dv = £. Given the vorticity 
2-form £, the Poisson equation A/3 = £ on the 2-forms is equivalent to Poisson equations for their 
respective components: A/3jj = Then /3 can be reconstructed component- wise by using the Green 
function: /%(g) = f Rn G(q,p) £y(p) p(p). 

For £ = £p and A/3 = £p = C ■ Sp, one has /3 VlV2 (q) = C J p G(q,p)pp(p) for the normal to P 
component of the potential /3, while other components are zero. Here pp is the volume form induced 
from 1" to P. 

By plugging this to the formula E(v) = \ J M £ A */3 and using *£p = C ■ pp we obtain 

C f C 2 f f 

E ( v ) = TT / Pvii* HP = / / G(q,p)pp(p)p P (q). 
1 Jp z JqeP JpeP 

The latter is a divergent integral, which can be regularized by considering 
E e {y) = 2 Jmn( v > u e) A 4 - Namely, given a point q € P replace the inner integral over P by the one 
over P e = {p £ P | ||g — p|| > e} by removing from P the e-neighborhood of q. Then one has 

E e(v) = %- I f G{q,p) p P (p) p P (q) . 

L Jq&PJp€Pc 
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As e — > the inner integral f &p G(q,p) p,p{p) increases as C n ■ lne, where the constant C n depends on 
dimension n only. Indeed, as p — >■ q one has G(q,p) = C n \\q — p\\ 2 ~ n ~ C n p 2 ~ n , where p is the geodesic 
distance from p to q in the membrane P. Then the integration of G(q, p) in the spherical coordinates 
over a small (n — 2)-dimensional punctured neighborhood ?7 €)a of radius a around the point q € P in 
the membrane P gives the integral 

f G(q,p)p P (p)~C n f p 2 ~ n p n ^dp = C n r p~ 1 dp = -C n lne + 0(l) 

as e — >■ 0. Then after the second integration over q G P the regularized energy E e (v) has the following 
asymptotics: 

E e (v)=C n -f Qne)np(q) + 0(1) = C n • In e • volume(P) + 0(1) as e^O, 
JqeP 

which completes the proof. □ 

Remark 6.5. One can see from the proof that the logarithmic singularity of the energy E e {v) comes 
from close points in P. To find the asymptotics one specifies a small parameter a giving the "range of 
interaction" and send e — > 0, while other pairs of points do not contribute to the leading term in the 
expansion of E t (v). This explains the term "localized induction approximation" (LIA). 

Rcnormalize time and regard H(P) := volume(P) as the new energy associated with fluid motions 
whose vorticity is supported on the membrane P. As we discussed above, this leads to the Hamiltonian 
dynamics of the membrane given by the skew-mean-curvature flow in any dimension. 

Note that in the regularization above one essentially uses only the symmetry and the order of 
singularity of the Green function G(q,p) as \\q — p\\ — » 0. The same asymptotics of the Green function 
holds for an arbitrary manifold M n and so does the energy regularization which results in H{P) = 
volume(P). 
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